Hierarchy in the Phase Space and Dark Matter Astronomy 
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We develop a theoretical framework for describing the hierarchical structure of the phase space 
of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the 
full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo 
model, the stable clustering hypothesis is applied for the first time here to the small-scale phase 
space structure. As an application, we show that the particle dark matter annihilation signal could 
be up to two orders of magnitude larger than that of the smooth halo within the Galactic virial 
radius. The local boost is inversely proportional to the smooth halo density, and thus is 0(1) within 
the solar radius, which could translate into interesting signatures for dark matter direct detection 
experiments: The temporal correlation of dark matter detection can change by a factor of 2 in the 
span of 10 years, while there will be significant correlations in the velocity space of dark matter 
particles. This can introduce 0(1) uncertainty in the direction of local dark matter wind, which was 
believed to be a benchmark of directional dark matter searches or the annual modulation signal. 



Among the favorite dark matter candidates are the 
weakly interacting massive particles (WIMP), yet to be 
detected in particle accelerators. Different experiments 
look for a signature of WIMPs directly as they pass 
through Earth and recoil off the atomic nuclei in lab- 
oratory detectors, or indirectly through the byproducts 
of their self-annihilation into standard model particles 
(such as photons, electron/positron pairs, or neutrinos). 
The cross-section for the (s-wave) self-annihilation is 
also fixed by the relic abundance: {cr ann v) ~ 10~ 31 — 
10 _26 cm 3 /s. Therefore, one would expect an astrophys- 
ical luminosity of: 
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where m x is the WIMP mass, Ei is the energy of the 
byproducts, dBi/dEi is the differential branching ratio 
into i (=7, e + e~ or vv) particles and <E> is fixed by the 
spatial density distribution of dark matter within the 
emitting region. 

All the factors in ([!]), except for $, are fixed by the 
particle physics model (e.g. jlj). Unfortunately the grav- 
itational potential of dark matter haloes says little about 
the contribution of small scale structure to $. The cold 
dark matter (CDM) primordial power spectrum predicts 
a large range of mass scales, from 10 12 — 1O 14 M down 
to 1CT 12 - 1CT 4 A/ Q @, far below the resolution limit of 
the present-day simulations at z = 0, which is at best 
10 4 M Q . In fact, simulations already see a nearly con- 
stant contribution to <!> per decade in sub-structure mass 
||, 9, suggesting a significant contribution from unre- 
solved structures. 

Extrapolation of the simulated properties of the sub- 
haloes to below the resolution limit, as well as assump- 



tions about the spatial/mass dependence of the boost 
due to sub-structure (e.g. f|, 0) could lead to signif- 
icant over/underestimations (especially for a non-scale- 
invariant linear power spectrum such as CDM). Similarly, 
the assumption of Maxwell-Boltzmann velocity distribu- 
tion in direct detection experiments of dark matter par- 
ticles could miss important phenomcnological signatures 
of clustering in the phase space. 

In this letter, we use the stable clustering hypothesis 
in order to predict the clustering of dark matter parti- 
cles in the phase space, and consider the implications for 
indirect and direct dark matter searches. The stable clus- 
tering hypothesis was first introduced by Davis & Peebles 
H as an analytic technique to study the galaxy correla- 
tion function in the deeply non-linear regime, and was 
subsequently applied into fitting formulae for non-linear 
correlation functions (e.g. (|, [Io| [n}). The hypothesis as- 
sumes that the number of neighbors within a fixed phys- 
ical separation becomes a constant (or pairwise velocity 
vanishes) on small scales, when the non-linear structure 
formation is completed. However, it became clear that 
this cannot be a good approximation on the scale of virial 
radius of CDM haloes (~ 1 Mpc today) as halo mergers 
and the subsequent tidal disruption can dissolve the old 
structures into newly formed haloes (e.g. [|l2|, |l3|). Nev- 
ertheless, the alternative analytic framework, known as 
the halo model, which has been extensively used in the 
literature over the past decade, misses the small scale 
structure of the haloes (i.e. sub-haloes, sub-sub-haloes, 
etc.) that do naturally form in hierarchical structure for- 
mation, and are ubiquitous in high resolution N-body 
simulations (e.g. 0]). 

In contrast, the merging and tidal activity should cease 
at scales much smaller than the virial radius of CDM 
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haloes, suggesting that the stable clustering regime might 
be achieved on small enough scales. In other words, while 
sub-haloes could lose a large fraction of their mass due 
to tidal heating/stripping, a small fraction could remain 
gravitationally bound (e.g. As the mcan/virial 

density of haloes drops as 1/time 2 , the gravitationally 
bound remnant will eventually become resilient to tidal 
disturbances. The stable clustering hypothesis can be 
trivially extended to the phase space where similar to 
real space it would predict that the number of particles 
within the physical velocity Av and physical distance Ar 
of a given particle does not change with time for small 
enough Av and Ar. 

In order to develop the stable clustering formalism in 
phase space, we start with the collisionless Boltzmann 
equation at the phase space coordinates, r + Ar, v + Av, 



df 



(r + Ar, v + Av,t) 



dt 

df , df 



df 



ill Or ( v + Av )-^-(V0 + VV0-Ar)=O.(2) 
The above equation in terms of the new function 

/ l (Ar,Av) = /(r l + Ar,v, + Av) (3) 
for particle i, can be re- written as 



df = 

dt 
where 
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■(VV0'Ar) =0 (4) 
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(5) 



The stable clustering hypothesis assumes that the above 
expression when averaged over the particles vanishes for 
small Ar and Av. If we assume that (/iVV0)p ~ 
(/i)p(VV^)p then a solution to (||), averaged over parti- 
cles is 

</) p = jjJ2h= F t Av2 + Ax, Ax fc (d,cW> P ] • (6) 



which is the most general time-independent solution with 
an isotropic velocity distribution, where F is an arbitrary 
function, and N is the number of particles in the phase- 
space volume of interest. We later show that averaging 
over particles and volume-averaging differ only by a con- 
stant. Next, we use the approximation that the potential 
is spherically symmetric and hence the above solution, 
using Poisson equation, can be re-written as 



</>p =Ks=F [(Av) 2 + 10(lff(£ s ) 2 (Ar) 2 ] 



(7) 



where £ s is the phase-space density at the formation time 
of the sub-structure, and we have used the fact that the 



post-collapse halo density is roughly ~ 200 times the crit- 
ical density at the formation time /i ~ 1 — 10% is 
the mean fraction of bound particles that can survive the 
tidal disruption period. 

In order to find the function F, we use the spherical 
collapse results 



is 



G 2 M(£ S ) ' 



(8) 



and also that the radius and velocity dispersion of haloes 



are related by (e.g. M): cr v 



10Hr v 



where H 



is the Hubble constant at the time of halo's collapse. 
Hence (JsJ) is roughly the phase space density of haloes 
that collapse at Hubble constant H(£ s ) and mass M(£ s ). 
The phase space volume of the collapsed halo, i.e. the 
volume of the constant-^ ellipsoid in (Q) is M/£ s and 
using (N) we have: 



[GM(Zs)Y 
10H(6) 



(9) 



Furthermore, the mass scale that collapses at a given 
cosmological epoch is characterized by: 



H 



-2/3 



v[M(£ 3 )] ~£ c ~ 1.7 



(10) 



where S c is the linear density threshold for the spherical 
collapse, while a[M] is the r.m.s. top-hat linear over- 
density at the mass scale M, and Hq is the present day 
Hubble constant. We remark that the effect of dark en- 
ergy will be a constant factor that could be absorbed in 
the definition of /x, since most sub-haloes have formed 
long before the era of dark energy dominance. 

Using the above results, the phase space correlation 
function takes the form: 

(/(ri,vi)/(r 2 ,v 2 )) 

~ — / d 3 r d 3 v /(r, v) /(r + Ar, v + Av) 

Jv 6 

i 

~ (/(r 1 ,v 1 ))(/(r 2) v 2 ))+M(/(r,v))e s (Ar,Av), 

(12) 

where we used the assumption of ergodicity to replace 
the ensemble average () by the volume average, in a 
given volume of the phase space Ve, while (f, v) are the 
mean values of (ri, vi) and (ra, v 2 ). The second term in 
Eq. ( |l2| ) dominates in the stable clustering regime, where 
|Av| = |vi - v 2 | < Aw tid and |Ar| = |r x - r 2 | < Ar tid , 
and Awtid and Ar t id characterize the tidal truncation ra- 
dius in the phase space. On the other hand, the first 
term dominates Eq. ([l2]) for large separations in the 
phase space, where particles are not correlated. In other 
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words, Eq. (H) is an interpolation between the stable 
clustering and the smooth halo regimes. 

The annihilation signal from gravitationally bound 
substructure is found by integrating (Q) within the 
phase-space stable clustering hypothesis: 

5$sub. y"d 3 xy"d 3 v 1 d 3 v 2 (/(x,v 1 )) M 6(0,v 1 -v 2 ) 

= /d 3 xp halo (x) /<J 3 Av^ 8 (0,Av) , (13) 



yielding the local boost factor of : 

n , s _ (^halo(x) 2 ) _ A* 
i>sub.(. x ) = 



(Phalo(x)) 2 
l-K 1 / 2 ^ ( Phalo(x 



9^3 V200/9 C 



Phalo(x) 
!/-M max 



1 d 3 Av.£ s (0,Av) 



j\/„ 



M- 2 d[M 2 a 3 (M)] (14) 



where we have used the stable clustering framework de- 
veloped in (0)-(|lO|) to substitute for £ s . 

The asymptotic form of CDM linear power spectrum 
on small scales is (e.g. |l7| ) P(k) cx fc" s_4 ln 2 (/c/fc eq ) , 
where n s is the primordial adiabatic scalar index, and 
fc oq ~ 0.56 fl m h 2 Mpc^ 1 is related to the comoving scale 
of the horizon at matter-radiation equality. Using this 
asymptotic form, we can find an analytic approximation 
for the boost factor: 



.B su b.(x) = /-« 



K 



( Phalo(x) 
V200pcrit,0 
f M cq 



£1 m ft 



0.27 x 0.7 



[( "eft, 8+3) 



In 1 / 2 



-K 



In 1 / 2 



cq 



M„ 



15) 



where M cq = (2.9 x 10 14 M )(f2 m /i 2 ) _2 , is the mass asso- 
ciated with the horizon scale at matter-radiation equal- 
ity, n e ff,8 — —1-68, is the logarithmic slope of the linear 
power spectrum at ~ 8/i _1 Mpc, and 



yy> 10 5 V 11 



„9 



y exp 



\{y 2 - i6)K 



(16) 



This analytic approximation (113) is within 40% of the 
exact integral in Q for 10^ 1 ^M Q < M < 10 10 Af o and 
0.95 < n s < 1.05 (where we compare with the fitting form 
of [pjj for the CDM power spectrum). For supersymmet- 
ric dark matter models, the minimum CDM halo mass 
from lO- 12 Af -lO- 4 M , yielding K ~ 10 5 



can range 
within a factor of 3. 

To calibrate the parameter (i, we can compare our re- 
sults to numerical simulations. Matching the expected 
boost per mass decade to that of "Via Lactea" mock 
Milky Way simulation which is around 0.1 for sub- 
halo masses 10 7 M Q < M sub < 1O 1O M , with their as- 
sumed cosmology, yields [i ~ 0.026 , which is comparable 
to the fraction of simulated halo mass found in (resolved) 
gravitationally bound sub-haloes (e.g. 5.3% in ||). 
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FIG. 1: The cumulative radial profile of DM annihilation: top 
dotted curve shows the contribution from the mean profile, 
while the lower dotted curves show the simulated contribution 
of sub-haloes more massive than 10 5 ' 6,7 ' 8 A'f© M. The solid 
curves show our analytic prediction, with the same cut-off's 
on Mmin, and a tidal cut-off at the high mass end. The dashed 
curve shows the cumulative mass profile for comparison. 



Fig.(Q) compares the prediction for the annihilation 
profile (solid curves) with numerical simulations of Q] 
(dotted curves), which shows reasonable agreement. To 
do this, we have cut off the integral (|l4|) at high masses 
by requiring that the formation density of bound sub- 
haloes must be lower than the local halo density. The 
crudeness of this criterion is most likely responsible for 
the sharp drop in the signal at small radii, relative to 
the numerical results. However, we note that (surviving) 
sub-solar mass subhaloes, that dominate the boost for re- 
alistic WIMP models, are much less affected by the tidal 
effects of the host halo (which was the motivation for the 
stable clustering approximation). Moreover, the simu- 
lated results do not include the effect of sub-sub-haloes, 
etc. which could be another source of difference with our 
prediction. 

Fig. (0) shows that not only the CDM hierarchy can 
significantly boost the total annihilation signal from the 
Milky Way (and dark matter haloes in general) , they can 
also affect the local density variance at the solar radius 
(B su h. ~ 2 — 6), which can have interesting implications 
for direct dark matter detection searches. We can quan- 
tify the latter through the temporal correlation of dark 
matter detection signal D(t), which traces the local den- 
sity of the dark matter halo at the solar system. For 
simplicity, we will assume that the solar system is mov- 
ing through the CDM hierarchy at v ~ \/2.5 x 250 km/s 
(assuming a singular isothermal mean phase space dis- 
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FIG. 2: Boost to the dark matter annihilation signal due 
to substructure with the formation mass larger than M m i n . 
The lower solid curve is the local boost factor at solar radius 
(assuming p = 10 pcrit), while the upper solid curve shows 
the mean boost estimated for the whole Milky Way halo (out 
to 16 times NFW scale radius H), where we have assumed 
n, = 0.96, a 8 = 0.82, Q m = 0.28, and h = 0.7. The dotted 
curve is our analytic approximation to the boost factor ([l5|). 
The dashed region shows the theoretical expectation for the 
mass cut-off I 



tribution). The two-point correlation function of D{t) 
simply measures the two-point correlation of CDM den- 
sity, by projecting ( |l2] ) into the real space, which is then 
modulated by the annual motion of earth around the sun. 
This can be done through a simple generalization of Ji4| ) 
to allow for finite separation in real space, and is shown 
in Fig. (||) for three different cut-offs of the CDM hierar- 
chy, assuming pcbm = lO'Vcrit within the solar system. 
We thus predict that, depending on the cut-off in the 
CDM hierarchy, the dark matter detection signal could 
gradually change by up to a factor of 2 within a ten year 
period. Potential measurement of this temporal correla- 
tion could shed light on the cut-off of the CDM hierarchy, 
which is directly related to the mass of CDM particles. 

To conclude this letter, we will turn to directional dark 
matter searches, which can directly probe the phase space 
of the CDM halo. These searches are now underway, and 
can potentially provide the first telescopes for the rich 
field of Dark Matter Astronomy. We can use our for- 
malism (Eqs. [7 TO) to predict the velocity space corre- 
lation function at the solar radius, which can be directly 
measured if CDM particles are detected in directional 
searches. This is shown in Fig. (|J). For simplicity, we 
assumed a singular isothermal sphere for the mean phase 
space density at the solar radius. The dimcnsionlcss cor- 
relation function is more prominent at higher velocities, 



FIG. 3: Temporal correlation function of the dark matter 
detection signal within the solar system (assuming pcdm = 
10 J p C rit) for three different cut-offs of the CDM hierarchy. 
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FIG. 4: DM phase space correlation due to bound subhaloes 
at solar radius. The three lines show velocities: 200, 400, and 
600 km/s in the Galactic frame. The shaded area shows the 
expected contribution from unbound substructure ]l^ ]. 



and shows a correlation length of 5-20 km/s in the ve- 
locity space. For comparison, the shaded area shows the 
expected level of correlation due to unbound substruc- 
ture, which is 10-30% at Av ~ 100 km/s pfj, and scales 



The fact that 



2-6 implies that the local CDM 



density may be dominated by small subhaloes with ran- 
dom velocities with respect to the Galaxy. This will 
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introduce an 0(1) uncertainty in the direction of local 
dark matter wind (or dipole), which was believed to be 
the benchmark of directional dark matter searches, or 
the annual modulation signal. Nevertheless, as we ar- 
gued above, the richness of structure introduced by the 
CDM hierarchy leads to novel observables/smoking guns 
for dark matter searches to aim for. 
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